Monte Carlo study of the transverse susceptibility in ordered arrays of magnetic 
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We present Monte Carlo simulations of the reversible transverse susceptibility (RTS) for a hexag- 
onal array of dipolar interacting magnetic nanoparticles with random anisotropy. RTS curves with 
the bias-field in-plane and out-of-plane are compared. With increasing temperature the RTS curves 
evolve from a three-peak {He, =LHic) structure to a double-peak and eventually a single-peak at the 
blocking temperature of the system. This trend is preserved for weak interactions. Dipolar inter- 
actions at low temperature are responsible for the suppression of the He peak in the out-of-plane 
geometry and its progressive merge to the Hk peak with decreasing interparticle separation in the 
in-plane geometry. The Hk peaks are located at higher field values in the out-of-plane geometry 
relative to the in-plane one. When the bias field lies in-plane (out-of-plane) the Hk peaks are shown 
to shift to lower (higher) field values with decreasing interparticle separation. The He peak shifts 
to lower field values in both geometries. Our results are compared with recent experimental findings 
in self- assembled arrays of Fe nanoparticles. 

PACS numbers: 75.50.Tt, 75.75.+a, 75.20.-g, 75.30.Gw 

Keywords: magnetic nanoparticles; self-assemble ; transverse susceptibility; dipolar interactions; Monte Carlo 
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I. INTRODUCTION 

Ordered arrays of single-domain magnetic nanopar- 
ticles (NPs) produced by synthetic methods and self- 
assembly have attracted a lot of research effort over the 
last decadei, motivated mainly by the wide range of po- 
tential technological applications that vary from high- 
density magnetic storage media 2 to high-sensitivity field 
sensors and logic devices^. Owing to the periodic ar- 
rangements of the NPs, their single-domain phase and 
their high monodispersity, they constitute ideal systems 
to gain basic understanding and possibly control of the 
role of interparticle magnetostatic interactions in their 
static and dynamical magnetic behavior. The issue of 
interparticle dipolar interactions has been addressed so 
far by a variety of experimental techniques including 
in most cases SQUID magnetometry and AC suscep- 
tibility measurements^, small-angle neutron scattering 
(SANS)^ and resonant magnetic X-ray scattering^ the 
last two methods being direct probes of magnetic cor- 
relations at the interparticle scale. These studies have 
provided ample evidence that the interplay between ran- 
dom anisotropy and dipolar interactions determine the 
magnetic behavior of the NP arrays. In particular ob- 
servations such as, different values between the in-plane 
and normal-to-plane remanence magnetization 7 , distri- 
bution of energy barriers with a larger width than the 
corresponding particle volume distribution-., flat field- 
cooled magnetization curves^, increase of the blocking 
temperature with increasing number of monolayersiSiii 
have been attributed to interparticle dipolar interactions. 
Model studies of the field and temperature dependent 
magnetization 7 ^ of NP arrays demonstrated the role of 
dipolar interactions and supported most of the experi- 
mental observations. 

In addition to the above mentioned techniques, the 



reversible transverse susceptibility technique was theo- 
retically introduced by Aharoni et ali& According to the 
predictions of Aharoni et ai, under conditions of coherent 
rotation of the magnetization (Stoner-Wolfarth model) 
the dependence of RTS on bias field exhibits three dis- 
tinct singularities, at the anisotropy fields ±Hk and at 
the coercive field H c , which render the method suitable 
to probe the magnetic anisotropy. Almost thirty years 
later the first successful realization of the methodic in 
BaFeO powders verified the theoretical predictions. The 
RTS technique is anticipated to be particularly suitable 
for analysis of the anisotropy of single-domain magnetic 
NPs, as the magnetization dynamics of those is satis- 
factorily described by the coherent rotation assumption. 
However, there are various factors that make difficult the 
identification of the RTS peak positions and the extrac- 
tion of the single-particle anisotropy strength. It has 
been previously demonstrated that particle size distribu- 
tion rounds all three peaksi^, orientational texture sup- 
presses the coercivity peaki^, interparticle interactions 
cause the coercive anisotropy peaks to merg o 16 : 17 , and 
thermal relaxation of the moments causes large shifts of 
the coercive peakiii^. Despite of these difficulties, the 
RTS technique was further developed and applied suc- 
cessfully to the analysis of the dynamical magnetic be- 
havior of magnetic NP assemblies These stud- 
ies have demonstrated the capability of the technique 
to probe efficiently the transition of the NPs from the 
blocked to the superparamagnetic regime by analyzing 
the field and temperature dependent RTS curves. Anal- 
ysis of the peak structure, provided evidence of the pres- 
ence of interparticle dipolar interactions in dense assem- 
blies of 7 - Fe 2 O^S. and CoSi^ NPs. The issue of in- 
terparticle dipolar interactions and the induced collec- 
tive magnetization dynamics was also studied in self- 
assembled arrays (SAA) of Fc 1 — and Fe-based 2 ^ NPs. 



2 



The authors demonstrated the anisotropy between in- 
plane and out-of-plane RTS and suggested the existence 
of an intermediate phase between the blocked and super- 
paramagnetic, where dipolar interactions dominate the 
dynamics of the assembly. 

Previous theoretical investigations of the RTS in dipo- 
lar interacting NP assemblies have been performed within 
the mean-field approximation for the local field and a ki- 
netic equation approach to the description of the thermal 
relaxation of the magnetic moments^ This approach 
although it predicted correctly the fast decay of coer- 
cive peak with temperature it failed to give the expected 
shift of the anisotropy peak with temperature. Micro- 
magnetic studies of the transverse susceptibility imple- 
menting the Landau-Lifshitz-Gilbert (LLG) equations of 
motions for dipolar coupled magnetic moments" have 
also appeared. The singularities of RTS in samples with 
uniaxial and cubic anisotropy were reproduced and the 
broadening of the peaks due to dipolar interactions was 
demonstrated. However thermal relaxation of the mo- 
ments and orientational randomness of easy axes was not 
considered in that work and therefore the results can- 
not compare directly to existing experiments on SAA 
of magnetic NPs. Thermal fluctuations effects on the 
transverse susceptibility have also been considered im- 
plementing the stochastic LLG equations to calculate the 
imaginary part of RTS. 25,26 Yuan and Victora 25 studied 
granular films, where exchange interactions dominate. 
On the other hand, Cimpoesu et al 26 studied rectan- 
gular arrays of dipolar coupled NPs with perpendicular 
anisotropy and identified the effects of dipolar interac- 
tions on the signal of the imaginary part of RTS. Despite 
the interesting conclusions demonstrated in the work of 
Cimpoesu et al regarding the interplay of size distribu- 
tion and dipolar effects, these can not be extended to the 
case of SAA of magnetic NPs, which are characterized by 
random anisotropy and a hexagonal arrangement of the 
NPs. 

The aim of the present work is to model the behavior 
of the field dependent RTS of ordered arrays of mag- 
netic NPs taking into account both thermal fluctuations 
and dipolar interaction effects. This is achieved imple- 
menting the Metropolis Monte Carlo (MC) algorithm. 
With this algorithm, the fluctuations in the the local 
field are treated exactly (up to the accuracy imposed by 
the calculation of the long-range dipolar forces) and the 
magnetization correlations, required to describe collec- 
tive magnetic behavior, are properly developed during 
the simulation. 27 Finally, the arrangement of the NPs on 
the triangular lattice is a crucial ingredient of the model, 
because the dipolar interactions have a well established 
ferromagnetic (FM) character in this geometryiLi 2 ^ In 
Section II we describe our model for the NP array and 
the simulation method for the calculation of RTS, in Sec- 
tion III we present numerical results for the in-plane and 
out-of plane RTS and finally in Section IV we discuss our 
results and make a connection to related experiments. 



II. MODEL AND SIMULATION METHOD 

We consider N identical spherical particles with diame- 
ter D forming a two-dimensional triangular lattice in the 
xy plane with lattice constant d, where d> D. The size 
dispersion of the NPs can be neglected to a good approx- 
imation, since in most self- assembled samples a very nar- 
row size distribution (a < 5%) is achieved^ The particles 
are assumed single domain, with uniaxial anisotropy in 
a random direction, and they interact via dipolar forces. 
The total energy of the system is 

E = 9 ±^ Rj- 

-k^Si-e^-h^Si-H) (1) 

where Si is the magnetic moment direction (spin) of par- 
ticle i, el is the easy-axis direction, and Rij is the center- 
to-center distance between particles i and j. Hats in 
Eq. indicate unit vectors. Three energy parameters 
enter Eq. QJ, namely (i) the dipolar energy g — m? /d 3 , 
where m — M S V is the particle moment, M s the sat- 
uration magnetization density and V the particle vol- 
ume, (ii) the anisotropy energy k = K{V , where K\ is 
the uniaxial anisotropy energy density and (iii) the Zee- 
man energy h = mH, where H is the applied dc field 
. The relative strength of the energy parameters en- 
tering Eq. l[TJl. the thermal energy t = ksT, and the 
treatment history of the sample determine the micro- 
magnetic configuration of the assembly. In all subse- 
quent results we scale all energy parameters by the sin- 
gle particle anisotropy energy (k = 1). The transition 
from single-particle to collective behavior is determined 
solely by the ratio of the dipolar to the anisotropy energy 
g/k = {■K/&){M 2 /K 1 ){D/df. The reported value* 7 -^ 
for fee or hep Co NPs are g/k = 0.2 - 0.4(D/d) 3 , while 
for the soft e-Co phase, higher values are expected^ For 
Fe NPs Farrell et al 4 report g/k = 1.54 {D/df and Pod- 
dar et aiiS report g/k — 2.8 (D/d) 3 . To compare with 
the experiments of Poddar et al we choose Fe NPs with 
D = 6.8nm at a separation d ~ 20nm, which is an esti- 
mate from their TEM images^ These values correspond 
to g/k = 0.1 in our simulations. 

Reversible transverse susceptibility measurements are 
performed with an ac measuring field (H ac ) perpendicu- 
lar to the dc bias field. The ac field is weak (~ 10 Oe) and 
its frequency lies in the rf regime (/ ~ 10 6 Hz)m&12£&. 
These experimental conditions allow the following ap- 
proximations to be adopted in the calculation of the RTS. 
First, since the amplitude of the measuring field is neg- 
ligible compared to the saturation field (~ 10 3 Oe), the 
calculation of the susceptibility is performed in the zero- 
field limit (H ac — 0). Second, the frequency dependence 
of the measuring field is neglected. This approximation 
is justified as long as the relaxation time of the NPs is 
large compared to the inverse frequency of the measur- 
ing field. Assuming that the NPs obey the Neel-Brown 
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FIG. 1: Sketch of the geometry used to calculate RTS. (a) In- 
plane geometry, used to obtain the parallel RTS (Xt)> an d (b) 
out-of-plane geometry, used to obtain the perpendicular RTS 
(Xt)- The NP superlattice is generated with basis vectors 
Si = (1,0) and a 2 = (5,^). 



model for thermal relaxation, their relaxation time is 
given by21 T = f~ l exp{-KV/k B T) with f ~ 10 w Hz 
and the blocking temperature = KV/ln(f a T m )kB- 
For RTS measurements the characteristic measuring time 
Tm = 2 • 10~ 5 s, and one obtains rf = 100 for 
T = 0.94T 6 and r/ = 10 for T = 1.16T&. Therefore, 
the static approximation is reasonably justified for tem- 
peratures up to ~ 20% above the blocking. In addition, 
our focus in the present work is on the static properties 
and their modification due to dipolar interactions, rather 
than on frequency dependent quantities, thus working in 
the static limit serves our purpose. For in-plane mea- 
surements we take the bias field H^c along the x-axis and 
for out-of-plane measurements along the z-axis. In both 
cases the measuring field is assumed along the y-axis (see 
Fig. [TJ. The values of the RTS (per spin) are calculated 
from the fluctuations of the magnetization, M y — ^2 Sf , 
as 
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where Xt ^ s the in-plane (out-of-plane) RTS (see 
Fie. [T)l. Simulations were performed for an ensemble of 
N=400 spins located in a simulation cell with dimen- 
sions L x = 20d and L y = 10v3d cut from a triangular 
lattice. For the interparticle interactions we used free 
boundaries along the z-axis and periodic boundaries in 
the xy plane that diminish undesirable in-plane demag- 
netizing effects arising from free poles. The dipolar inter- 
actions were summed to infinite order in-plane, using the 
Ewald summation method for a quasi-two-dimensional 
system; 30 . For the simulation of the magnetic configura- 
tion under an applied field H and at finite temperature 
T we used the standard Metropolis MC algorithm with 
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FIG. 2: Temperature dependence of remanence (left) and co- 
ercivity (right). Circles: non-interacting particles. Squares: 
interacting (g = 0.1) particles and in-plane field. Triangles: 
interacting (g = 0.1) particles and out-of-plane field. 



single-spin moves^i According to this algorithm, a spin 
is chosen at random and it is rotated by a small an- 
gle. This is achieved by varying the Cartesian coordi- 
nates (Sf, Sf , Sf) of the moment randomly in the inter- 
val (—(5,(5), where < 5 < 1, and renormalizing its new 
magnitude to unity. The new configuration is accepted 
with probability equal to mm[l, exp{— AE/ksT)], where 
AE is the change of the total energy. The value of S is 
adjusted such that approximately 50% of the attempted 
moves are accepted. Starting from a chosen spin config- 
uration, the initial 10 3 MC steps per spin (MCS) were 
used for relaxation of the system towards equilibrium 
and thermal averages were calculated over the subsequent 
10 4 MCS, allowing 10 MCS between sampling events to 
achieve statistical independence. The results were aver- 
aged over N c = 100 samples with different realizations 
of the random axes distribution and the thermal fluctu- 
ations. In all results presented below a sweep of the bias 
field from negative to positive values is performed with a 
step of Ah — (l/30)fc. As a test of convergence a sweep 
in the opposite direction was performed in certain cases, 
but the deviations obtained in the positions of the peaks 
of RTS were well within the statistical errors. 



III. 



NUMERICAL RESULTS 



A. Remanence and coercivity 

We discuss first the characteristics of the hysteresis 
loop for a NP array. Dipolar interactions on a trian- 
gular lattice favor a ferromagnetic ground state and in- 
troduce an in-plane anisotropyiiSi 3 ^ 3 " The interaction- 
induced anisotropy produces a different magnetic behav- 
ior for in-plane and out-of-plane directions of the bias 



field. In Fig. we show the temperature dependence of 
the remanence (m r — M r /M s ) and coercivity (h c ) of a 
NP array with and without dipolar interactions. We no- 
tice that the in-plane remanence (m|) is enhanced while 
the out-of-plane remanence (m^r) is suppressed relative 
to the values for the non-interacting array. Further- 
more, the in- plane coercivity (/z|) is clearly enhanced 
due to the extra barrier provided by the dipolar inter- 
actions, while the out-of-plane coercivity (hj:) shows a 
weaker dependence on the interaction strength. For the 
non-interacting assembly the remanence and the coerciv- 
ity exhibit vanishingly small values, m r (t)/m r (0.05) < 
1% and h c (t)/h c {0 -05) < 0.5%, at temperatures above 
t ~ 0.14, which is therefore the blocking temperature 
for the non-interacting assembly (t°)£^ On the other 
hand, for the interacting assembly similarly small val- 
ues, m r (t)/m r (0.05) < 2% and h c (t)/h c (0.05) < 0.5%, 
are obtained above % ~ 0.17. At temperatures above 
the remanence and the coercivity are enhanced due 
to dipolar interactions. This result defines an interesting 
temperature regime, t® < t < tb, in which the thermal en- 
ergy overwhelms the random anisotropy barrier and the 
dynamics of the NPs is governed by dipolar interactions. 
The latter have a magnetizing effect at temperatures be- 
low i ,, as evidenced by the enhancement of mJL We refer 
to temperatures in the range t® < t < tb as the super- 
ferromagnetic (SFM) regime. Notice also that the data 
presented in Fig. 2 indicate that the SFM regime, is se- 
riously suppressed in the out-of-plane geometry, due to 
the strong demagnetizing character of the dipolar inter- 
actions in this geometry. Further down in this work we 
will discuss the behavior of the RTS in the SFM regime. 



B. Thermal evolution of transverse susceptibility 
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The field-dependent RTS of a non-interacting assem- 
bly is shown in Fig. |3| for various temperatures that ex- 
tend from low (t — 0.05) to high values (t = 0.50). 
At low temperature, the RTS curve shows two broad 
peaks at symmetric positions that correspond to the 
anisotropy field ±hx and a third sharper peak at the 
coercive field h c . The presence of three peaks in the low- 
temperature RTS curve is in accordance with the theo- 
retical predictions of Aharoni et ai— and previous nu- 
merical studiesii£*i£ Notice however that the anisotropy 
peaks appear at hx ~ ±1.0 which is lower than the 
zero-temperature value hx = ±2, because of the non- 
zero temperature at which the calculation is performed. 
Indeed, for an isolated NP, the anisotropy peak is cen- 
tered at the field required for an irreversible switch of 
the particle's moment. Therefore, the thermal energy 
assists the switch and causes a shift of the anisotropy 
peak towards lower field values. In previous numerical 
calculations of RTS at finite temperature^, the authors 
implemented a rate equation approach and found that 
the positions of the anisotropy peaks are insensitive to 
the temperature. This result does not agree with RTS 



FIG. 3: Evolution of RTS with temperature for a non- 
interacting (g — 0) array. The magnetic field is swept from 
negative to positive values. 



measurements in NP assemblies that show a clear down- 
shift of the anisotropy peaks with temperatur e) 19 ! 20 . Our 
approach, that goes beyond the mean field description 
of the thermal fluctuations of the magnetization predicts 
the expected downshift of the anisotropy peaks with tem- 
perature. As shown in Fig. [21 the coercivity decays fast 
with temperature and in Fig. |3| we show that already 
for temperatures t = 0.10 — 0.11, which lie well below 
the blocking (i° = 0.14) the coercivity peak of RTS is 
suppressed, leading to a pair of asymmetric anisotropy 
peaks. The broader among the two anisotropy peaks 
appears at a positive field, namely at a strong enough 
field to cause reversal of the magnetization. The forma- 
tion of a broad anisotropy peak due to its merge with 
the coercivity peak has also been observed in RTS mea- 
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surements in Fe-based2£ and Co—i^ NP assemblies. The 
fast downshift of the h c peak relative to the slow shift of 
the Hk peaks was also found in previous numerical stud- 
ies of RTS and suggests that extraction of the h c field 
from RTS measurements is much more difficult than of 
the }ik field, even for monodisperse and extremely dilute 
samples^ii^ When the temperature is increased above 
t = 0.11 (Fig|3}, the anisotropy peaks rise and merge 
to a single narrow peak that subsequently gets broader 
and lower above t = 0.14. The susceptibility at zero field 
(Fig. |3J) exhibits its maximum value at t = 0.14, which 
is the blocking temperature of the assembly. Above the 
blocking , the zero-field RTS drops with temperature ac- 
cording to the Curie law (x ~ l/^ 1 )- In summary the 
thermal evolution of the RTS curves, shown in Fig. |3 
is characterized by a three-peak mode at very low tem- 
perature followed by to a double-peak mode due to the 
merge of the coercivity and anisotropy peaks at interme- 
diate temperature, and finally, a single-peak mode above 
the blocking temperature (t > t®). 

When weak interparticle interactions (g — 0.1) are 
switched on, they do not modify the overall trend of 
the RTS curves with temperature, namely the transition 
from three-peaks to a double-peak and finally to a single- 
peak (Fig. |2J| . However, even weak interactions introduce 
noticeable differences between the in-plane susceptibility 

(Xt) an( i the out-of-plane one (Xt)- Starting from the 
low-temperature regime (f = 0.05) we notice that the 
h c peak is suppressed in the Xt curve while it is quite 
pronounced in the Xt curve - The in-plane anisotropy 
induced by dipolar interactions transforms the spheri- 
cal distribution of the easy axes directions to a quasi 
two-dimensional distribution that causes suppression of 
the hjr peak. A similar suppression of the h c peak was 
demonstrated by Hoare et at— due to easy axes texture 
in a non-interacting assembly. For dipolar coupled ar- 
rays, the interparticle interactions produce a similar ef- 
fect to texturing via the long range demagnetizing field 
that forces the moments to stay in-plane. At all tem- 
peratures the Xt curve is broader than the Xt curve i 
because of the larger saturation field in the out-of-plane 
geometry. The slow saturation arises due to the in-plane 
anisotropy induced by dipolar interactions which renders 
magnetically harder the z-axis relative to the x,y-axes. 
A further difference occurring between the two measur- 
ing geometries at low temperature is that the anisotropy 
peaks are located at higher fields in the case of out-of- 
plane geometry {hj^ > h^ K ). This feature is in agreement 
with experimental observations in Fe NP arrays^ 

A physical argument to justify this behavior is the fol- 
lowing : The position of the anisotropy peak corresponds 
the field required for an irreversible switch of the parti- 
cle's moment over the uniaxial anisotropy barrier. Dipo- 
lar interactions induce an easy-plane anisotropy addition- 
ally to the uniaxial one. Therefore, the in-plane switching 
is facilitated by the presence of more than one easy di- 
rections and consequently the switching field is reduced, 
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FIG. 4: Evolution of RTS with temperature for a dipolar 
interacting (g — 0.1) array. The magnetic field is swept from 
negative to positive values. Closed circles : in-plane field. 
Open circles : out-of-plane field. 



leading to a downshift of the corresponding RTS peak. 
On the contrary, out-of-plane switching is obstructed by 
the additional easy-plane that tends to keep the moment 
inside the xy-plane and the out-of-plane switching field 
is enhanced. Eventually, as the temperature rises ther- 
mal fluctuations overwhelm the dipolar anisotropy and 

the difference between the x\ anc ^ Xt curves is reduced, 
as shown in Fig. 0] for t = 0.30. 

To extract the blocking temperature from the suscepti- 
bility curves, we plot in Fig.lSJthe zero- field susceptibility 
Xt(H = 0) as a function of temperature. Provided that 
the peak of xt(T;H — 0) occurs at the blocking tem- 
perature of the system, we obtain from Fig. [3] the values 
t® ~ 0.14 and ff, ~ 0.17. These values are in satisfac- 
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FIG. 5: Temperature dependence of the zero-bias RTS. Cir- 
cles: g = 0. Squares: g = 0.1 and in-plane geometry. Trian- 
gles: g = 0.1 and out-of-plane geometry. 
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FIG. 6: Field dependent RTS in the superferromagnetic 
regime (t = 0.15). Open circles : g = 0. Closed circles : 
g = 0.1 and in-plane geometry. Field step : (a) Ah — (l/30)fc 
and (b) Ah = (l/300)fc. 
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tory agreement with the values obtained earlier from the 
temperature dependence of the remanence and coerciv- 
ity (Fig. |2J) . Notice that in contrast to the experimental 
situation, where magnetization and RTS measurements 
have different measuring times {tsquid ~ 10 2 s, trts ~ 
10~ 5 s) 10 ' 22 leading to different estimates of the blocking 
temperature, in our case both magnetization and suscep- 
tibility data are obtained during the same Monte Carlo 
relaxation process and therefore correspond to measure- 
ments at the same time scale. 

A point concerning the thermal evolution of the RTS 
curves across the blocked to superparamagnetic regime 
is next. The issue is whether the transition from the 
two-peak to the single peak mode occurs below or at 
the blocking temperature. Previous calculations for non- 
interacting assemblies^ point to the second scenario, 
namely that the double peak structure is preserved right 
up to the blocking temperature. However, in recent ex- 
periments in densely packed Fe NP arrays Poddar el ati^ 
argued that the anisotropy peaks merge at a temperature 
(Tcross) that lies below the blocking temperature of the 



FIG. 7: Low-temperature (t = 0.05) RTS for variable dipolar 
strength. Closed circles : in-plane geometry. Open circles : 
out-of-plane geometry. The arrows in the out-of-plane data 
indicate the position of the coercive held. 



array. They supported this observation by arguing that 
the merge of the peaks indicates the overcome of the uni- 
axial anisotropy barrier by the thermal energy, while the 
maximization of the RTS peak corresponds to the over- 
come of the interaction-induced barrier by the thermal 
energy. 

We examined this argument by performing more re- 
fined simulations in the field range around the peak. In 
particular, we have performed a high-resolution simula- 
tion at a temperature t — 0.15, that lies in the SFM 
regime (see Fig. [SJl, using an order of magnitude smaller 
field step dh = (l/300)fc. The choice of temperature is a 
typical value in the SFM regime, in which, as discussed 
earlier, thermal fluctuations are adequate to overcome 
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irreversible switching of the moments. On the contrary, 
in the out-of-plane geometry the anisotropy peak shifts to 
higher fields with increasing coupling strength since dipo- 
lar interactions inhibit the irreversible switching along 
the z-axis. The outcome of the interactions for the or- 
dered array is to produce higher values for the out-of- 
plane anisotropy field than the in-plane one {h-t > h)l), 
which is in agreement with the RTS measurementsi&n 
arrays of Fe NPs. 



FIG. 8: Variation of the (a) coercivity peak and (b) anisotropy 
peak position with the dipolar strength at low temperature 

(t = 0.05). Closed circles : in-plane geometry. Open circles : jy CONCLUSIONS 

out-of-plane geometry. 



the uniaxial anisotropy but the system remains ferro- 
magnetic due to dipolar interactions. The results are 
shown in Fig. EJ where it is demonstrated that the dou- 
ble peak structure for the interacting array is preserved at 
this temperature, while the non-interacting array shows 
a perfectly broad single peak. Repeating the simulations 
at higher temperatures but still below the blocking tem- 
perature (t < 0.17) we reached a double-peak curve with 
a closer distance between the peaks and reduced asym- 
metry. Therefore our simulations show that in the SFM 
regime (t° < t < tt) the double peak structure of the RTS 
curve is preserved and it transforms to a single peak at 
the blocking temperature, similarly to the case of non- 
interacting assemblies. 



C. Evolution of transverse susceptibility with 
particle spacing 

Varying the interparticle spacing (d), modifies the 
strength of the dipolar coupling as g ~ 1/d 3 . We show 
in Fig. [7] the in-plane and out-of-plane susceptibilities 
at low temperature (t — 0.05) for increasing dipolar 
strength. The most striking difference between the two 
geometries is that as the interaction strength increases, 
the anisotropy peaks of x\ shift to lower fields while those 
of Xt shift to higher values. The peak positions for the 
two geometries are plotted in Fig |S| as a function of the 
dipolar strength. 

The dependence of the anisotropy peaks on the cou- 
pling strength has been discussed above and is attributed 
to the easy-plane anisotropy induced by the dipolar in- 
teractions. A linear decrease of the anisotropy field with 
packing density has been previously reported in ran- 
domly packed magnetic NP»5^. Our simulations indi- 
cate that the same behavior is expected in ordered ar- 
rays and the in-plane geometry. The similarity between 
these two situations lies in the fact that in both cases the 
interaction-induced anisotropy reduces the barrier for an 



We studied the reversible transverse susceptibility in 
hexagonal arrays of anisotropic magnetic nanoparticles 
using Monte Carlo simulations, in order to investigate 
the effects of interparticle dipolar interactions. We found 
that below the blocking temperature dipolar interactions 
are responsible for a series of characteristic features in 
the RTS curves. Namely, (i) the suppression of the He 
peak of Xti ( u ) the location of ±Hk peak of Xt a ^ 
higher fields than the corresponding peaks of x\i (iii) t ne 
downshift (upshift) of the H'L (H^) peak with increasing 
dipolar strength or equivalently, decreasing interparticle 
spacing, (iv) the slower saturation with bias field of Xt 
relative to x\- These results are in agreement with re- 
cent measurements of in-plane and out-of-plane RTS in 
ordered arrays of Fe nanoparticles With respect to the 
thermal evolution of the RTS curves, we showed that 
with increasing temperature both the Hk and He peaks 
shift to lower field values and merge to a single peak 
leading to a double peak structure. Dipolar interactions 
are shown (Fig. to make the H c peak more persis- 
tent to thermal fluctuations. At the blocking tempera- 
ture of the system the anisotropy peaks merge and the 
zero-field RTS is maximized. The same thermal evolu- 
tion of the peaks is followed in weakly interacting arrays 
and agrees with the experimental observations in Fe NP 
arrays*^ Therefore, our simulations support the capabil- 
ity of RTS measurements to provide useful information 
related to the dynamical state of interacting nanoparticle 
assemblies. A minor difference between our simulations 
and measurements on SAA of Fe nanoparticles ^ is that 
experimentally the transition from a double-peak to a 
single-peak structure of RTS was observed at a tempera- 
ture T cross < Tf, and attributed to a transition to a super- 
ferromagnetic regime. Our simulations showed that the 
double peak structure persists up to T&. Possible reasons 
for this discrepancy should be sought for in the limita- 
tions imposed by the instrumental resolution or in the 
static approximations adopted in our model. This point 
requires further investigation. 
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